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Abstract 

In the context of modified theory of gravity (/(i?) gravity) we try 
to study the conditions needed for vaUdity of the generaUzed second 
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qh; 1 Introduction 

{3JQ| Recently, motivated by astrophysical data which indicate that the expansion 

of the universe is accelerating [1], the modified theory of gravity (or f{R) 
^SJ ■ gravity) which can explain the present acceleration without introducing dark 

I energy, has received intense attention [2]. 

Determining thermodynamic parameters of an (accelerated) expanding 
. universe and verification of the first and the second law for different cosmo- 

I logical horizons [3] ; investigating the relation between dynamics and thermo- 

0^ ' dynamics of the universe [Ij ; studying the conditions required for validity of 

the generalized second law (GSL) [5J, and so on, have also been the subjects 
i of many researches in recent years. 

In the modified theory of gravity, instead of Friedman equations we must 
utilize modified Friedmann equations which may include the powers of Ricci 
scalar as well as its time derivatives. Besides, the relation of the entropy to 
I the area of the horizon is also different with Einstein theory of gravity. So it 

is of interest to see how, in the framework of f{R) gravity, thermodynamic 
properties of the universe may be modified. 

In this paper we try to find necessary conditions for validity of GSL in 
the framework of f{R) gravity in Friedmann Robertson Walker (FRW) uni- 
verse. The future event horizon is taken as the horizon of the universe and 
the temperature is supposed to be proportional to the Gibbons-Hawking 
temperature [6]. The range of the proportionality constant will be deter- 
mined through some examples. Using the relation between the entropy as- 
signed to the horizon and its area derived from the Noether charge method. 
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we study the behavior of horizon entropy with respect to the comoving time. 
We assume that the horizon is in thermal equihbrium with its environment 
which is filled with perfect fluids. Then from modified Friedmann equations 
and the first law of thermodynamics, the time derivative of the fluid entropy 
and subsequently the time derivative of the total entropy is determined. At 
the end we elucidate our results via two examples. 
We use the units h = c = G = kB = l- 



2 Thermodynamics and GSL in FRW universe in 
modified gravity 

The action of modified theory of gravity with the inclusion of matter is given 

by 

S=r^j d^x^fiR) + S^, (1) 

where Sm is the matter action, R is the Ricci scalar curvature and f{R) is 
an arbitrary real function. Variation of the action with respect to the metric 
gives 

Rf.uf{R) - lgf.uf{R) + g^uafiR) - v^v,/' = 8^r-, (2) 

where the prime denotes the derivative with respect to i?, and TJJl is the 
energy-momentum tensor of the matter fields. For spatially flat FRW metric 
with scale factor a{t): 

ds^ = -dt^ + a^{t){dx^ + dy'^ + dz^), (3) 

Eq.@ yields 

Sttp = iM-3{H + H'-H^^)fiR) 

SvrP = -l^ + {H + 3H'-^^-2H^^)f'{R). (4) 

The Hubble parameter is given hy H = d/a, and the Ricci scalar is obtained 
as i? = GH + 12H^. The over dot indicates the derivative with respect to 
the comoving time t. p and P are the density and the pressure of the matter 
which behaves clS Si perfect fluid at large scale: 

T;:i = {p + P)U^U, + Pg^,, (5) 

where is the four velocity of the fluid. 

The radius of the future event horizon, Rh, is given by 
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Note that the event horizon exists, Rhit) G when the above integral 
converges, i.e., dt' / a{t') £ 3f?+. If at a time denoted by ts, the Big Rip 
singularity occurs, we must replace oo by tg in the integration. Using Eq.([6|), 
we can verify that satisfies the fohowing equation 

Rh = HRh - 1. (7) 

In the context of Einstein theory of gravity, the entropy of a black hole 
is given by Bekenstein-Hawking relation [7] 

Sh = ^, (8) 

where A is the area of the event horizon. In the same way one can assign an 
entropy to the cosmological future event horizon whose area is A^, = 47ri?^. 
This entropy which is given hy Sh = may be regarded as a measure of 
information hidden behind the horizon. In de Sitter space-time Rh = -jj 
and the future event horizon becomes the same as the de Sitter (Hubble) 
horizon. In this space time the temperature, which is dubbed as Gibbon 
Hawking temperature [6j, can be determined in terms of horizon radius as 

rp _ H_ 

~ 2n- 

In f{R) gravity, Noether charge method can be used to obtain the hori- 
zon entropy [8j 
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Sh = -Jj'iR)dA, (9) 

where the integration is taken over the surface of the horizon, A. So in 
FRW universe, where the scalar curvature is spatially constant, the entropy 
is obtained as 

5. = ^. (10) 

In the following we choose the future event horizon as the horizon of the 
universe. Differentiating Eq. (jlOp with respect to the comoving time gives 

Sh = 27rRhRhF + vriJ^F. (11) 

We have defined F = f'{R). In Einstein theory of gravity the above equation 
reduces to Sh = 2TTRhRh- Note that in a super-accelerated universe defined 
by if > 0, Rh is decreasing: Rh < [9], therefore Sh < 0. But in modified 
theory this is not the case and, depending on the function f{R), one may 
have Sh > 0. Consider the model f{R) = aR"^, a,m G For H > we 
have -R > 0, hence Sh >0 leads to 

It is clear that for a = m = 1 the above inequality cannot be satisfied in a 
super-accelerated universe. 
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If Sin is the entropy of the matter inside the horizon, then the first law 
of thermodynamics states 

TdSin = dE + PdV = {P + p)dV + Vdp. (13) 

By taking V = ^irRf^, we arrive at 

TS'in = MP + p)RlRh + ^T^RhP. (14) 

From Eq.([7|), and energy conservation relation 

p + miP + p) = Q, (15) 

we can write Eq. (jl4p in the form 

TS^n = -4vr(P + p)Rl. (16) 

Sin is a decreasing (increasing) function of time when w > {<) — 1, where 
w = ^ \s the effective equation of state (EOS) parameter of the perfect fluid 

filling the universe. In the model f{R) = R, Eq. (fT6|) becomes TSin = HRf^, 
hence Sin > is satisfied when H > 0. Note that in f{R) theory of gravity 
we have 

-4HF + 2HF -2F 
w = -l + ? ^ . 17 

f{R)-Q{H + m-Hj^)F 

Hence in contrast to Einstein theory of gravity, in f{R) models, H > does 
not requires w < —1 and H > is not a necessary condition for Sin > 0. 

Temperature of the horizon, T, which is taken the same as the fluid 
temperature, is independent of the gravity theory that leads to the horizon 
geometry (see e.g. Ref. (TOj)- In fiR) gravity, like Einstein theory of gravity, 
and in the absence of a well defined temperature for cosmological horizon, 
we assume that T is proportional to Gibbons-Hawking temperature 

-If- 

6 is a real constant. In a de Sitter space we must take 6 = 1. Indeed 
the parameter b shows the deviation from Gibbons-Hawking temperature. 
Inserting the modified Einstein- Friedmann equations (Eqs. (jl])), and (jlSp . 
into Eq. results in 

The total entropy of the universe, denoted by S, is the sum of the matter 
entropy inside the horizon. Sin, and Sh- S = Sin + S^. The generalized 
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second law (GSL) states that the total entropy is not a decreasing function 
of time: 

S = Sin + Sh>0. (20) 
Following our previous results, this leads to 

5 = .Hj((l-i)F + 2(| + ^)F+i^F)>0, (21) 

which is more complicated with respect to the model f{R) = R in which 
Eq. (j2ip reduces to the following simple inequality 



(HbRh) > 0. (22) 

To illustrate our results let us consider some examples. As the first 
example consider a quasi de Sitter FRW space time, defined by 

H = Ho + Hoet + 0{e^), e := « 1, e = 0{e^). (23) 
In this space time the future event horizon is given by 

RH = j^{l-§^) + 0{e'), (24) 

leading to 

Rh ^ -§^. (25) 
Therefore ^ — a-nd S > reduces to 

S ^ T^Rl (l -1){F- ^F) > 0. (26) 

Consider the model 

f = I3R + aR"", /3,a,me ^. (27) 

By using R ~ 2AHH and R ~ 24ij^ , one can verify that the GSL is satisfied 
up to the order ©(e^), provided that 

(1 - ^)H(am{m - 2)(12i/2)™-i - /?) > 0. (28) 

For a = and /? = 1, corresponding to Einstein theory of gravity, the above 
equation reduces to 

-(l-^)^gO. (29) 

showing that for a super-accelerated universe {H > 0), ^ is satisfied if 
6^1 and for quintessence phase {H < 0), GSL is respected when 6^1- 
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If /5 = and a > 0, S* ^ is satisfied when am{m — 2)(1 — \)H ^ 0. So 
if > and m > 2, we must have & ^ 1- Therefore in modified gravity in 
contrast to the Einstein theory of gravity, a super-accelerated FRW model, 
which departs slightly from de Sitter space, may have temperature greater 
than Gibbons-Hawking temperature and meanwhile respects the GSL. 

As another example consider again the f{R) gravity model defined by 
Eq. ([27]). Assume [ll] 

a = ao{ts-ty,n>Q. (30) 

In this model 

rr_ ^ _ \2v? + 6n 

Hence Eq. ([50]) describes a super accelerated FRW universe, H > 0, with a 
Big Rip singularity at t = ts- The future event horizon radius is 

Rk = '-^. (32) 
n + 1 

To show that Eq. ()30p may be a solution of Eq. ^ in /(i?) gravity, 
one can assume that the perfect fluid inside the event horizon includes two 
components with EOS: Pi = 'yipi and P2 = 72^2, where the constant 7's 
are the EOS parameters. These components satisfy 

Svrpi = ^-3P{H + H^) 

SttPi = -^+p(H + 3H^), (33) 



and 



'2i dt 



rv ^ rl^ rl 

8ttP2 = -^ + ma{H + 3H^ -^-2H^)R"'-\ (34) 
2 dt'^ at 

Besides, each component {i = 1, 2) satisfies the energy conservation equation 

Pi + 3H{-ii + l)pi = 0. (35) 

Using Eqs. (1330 . (I34p and Eq. ()35p and after some calculations it may 
be verified that Eq. @ is satisfied in the case (|27p provided that: 

^ = ^ = (36) 

which results in 71 < 0, and 72 < (>) — 1 if m > (<)0. Using Eq. ([2T]) . one 
can verify that in this model GSL is respected only for times satisfying 

-ma{{rn - 2){h - l)n + 2m^ + 1 - 3m) (^6n(l + 2n))'" ^ x 

(t, - t)2'"-2 + /3n(6 - 1) < 0. (37) 
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In Einstein theory of gravity {a = 0, (3 = 1), S'>Ois satisfied only when 
b < 1 , in this case 5 > holds \/t < ts- In the modified theory, depending on 
the values of m, n, a, and (3, S > holds only for times belonging to special 
domain specified by the Eq. (j37p and we may have 5 > while b > 1. We 
can also restrict the values of the parameters to specific domains such that 
GSL holds for Vt < ts, e.g., for ma [{m - 2)(b - l)n + 2m^ + 1 - 3m] > 0, 
GSL holds Vt < ts provided that - 1) < 0. 

It seems that GSL does not hold near the Big Rip singularity, t ~ t^, 
for P{b — 1) > 0. This may be related to the fact that in our classical com- 
putation we have ignored the contribution of the radiation energy density, 
generated by semiclassical particle creation from Rindler horizon near the 
Big Rip [T2], in the total entropy. 

For an adiabatic expansion (S = 0), the following equations hold: 

-ma((m - 2)(6 - l)n + 2m'^ + 1 - 3m) (6n{l + 2n)^'^ ^ = 0, 
Pn{b - 1) = 0, (38) 

which in Einstein theory of gravity (/3 = 1, a = 0), infers 6=1. Therefore, 
in this theory, during an adiabatic expansion the temperature is the same 
as Gibbons-Hawking temperature. While in the modified theory of gravity 
(a 7^ 0), in order to have S" = 0, (m — 2)(6 — l)n -|- 2m^ + 1 — 3m = and 
P{b — 1) = must be satisfied. Thereby /? 7^ leads to 6 = 1 and m = 1/2; 
and (3 = implies b = 1 + ^r^^^^^r^ ! so depending on the values of m 

Til { Tic \ 

and n, in adiabatic expansion, the temperature may be less or more than 
the Gibbons-Hawking temperature. 
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